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Abstract 

We reevaluate the hadronic contribution to the running of the QED fine structure constant a{s) 
ats = M|. We use data from e~^e annihilation and r decays at low energy and at the qq thresh- 
olds, where resonances occur. Using so-called spectral moments and the Operator Product Ex- 
pansion (OPE), it is shown that a reliable theoretical prediction of the hadronic production rate 
R{s) is available at relatively low energies. Its application improves significantly the precision on 
the hadronic vacuum polarization contribution. We obtain Aahad(-^z)~ (277.8 it 2.6) x 10~^ 
yielding a~i(M|) = 128.923 ± 0.036. Inserting this value in a global electroweak fit using 
current experimental input, we constrain the mass of the Standard Model Higgs boson to be 
MHiggs = 129+^^^ GeV. Analogously, we improve the precision of the hadronic contribution to 
the anomalous magnetic moment of the muon for which we obtain a^"^ = (695.1 ± 7.5) x 10-^°. 
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1 Introduction 



The running of the QED fine structure constant a{s) and the anomalous magnetic moment 
of the muon are famous observables whose theoretical precisions are limited by second 
order loop effects from hadronic vacuum polarization. Both magnitudes are related via 
dispersion relations to the hadronic production rate in e~^e~ annihilation, 

D/ \ (Ttot(e+e--> hadrons) 3s / + - uj n 

R{s) = — — ^ = rdtotie^e hadrons) . (1) 

While far from quark thresholds and at sufficiently high energy ^/s, R{s) can be predicted 
by perturbative QCD, theory may fail when resonances occur, i.e., local quark-hadron 
duality is broken. Fortunately, one can circumvent this drawback by using e~^e~ annihi- 
lation data of -R(s) and, as recently proposed in Ref. hadronic r decays benefitting 
from the largely conserved vector current (CVC). 

There is a strong interest in the electroweak phenomenology to reduce the uncertainty 
in a(M|) which at present is a serious limit to further progress in the determination of 
the Higgs mass from radiative corrections in the Standard Model. The most constraining 
observable so far has been sin^Qw obtained from leptonic asymmetries at the Z pole with 
an achieved precision of (Asin^Ovv)cxp = 0.00022 0. The uncertainty from the currently 
used a(M|) value translates into (Asin^Gvi/)^ = 0.00023 0, justifying the present work. 

In this paper, we extend the use of the theoretical QCD prediction of R{s) to energies 
down to 1.8 GeV. The reliability of this approach is justified by applying the Wilson 
Operator Product Expansion (OPE) 0] (also called SVZ approach [Q) and fitting the 
dominant nonperturbative power terms directly to the data by means of spectral mo- 
ments. An analogous approach has been successfully applied to the theoretical prediction 
of the T hadronic width, Rr, at Mr ~ 1.8 GeV in order to measure the strong coupling 
constant ^.(M^) [g ^ |, |. 

The paper is organized as follows: first a brief overview over the formulae used is 
given, then the spectral moments are defined and evaluated experimentally, and the cor- 
responding theoretical evaluation is fitted. On this basis, theoretical predictions of the 
hadronic vacuum polarization contribution to the running of a(M|) and to the anomalous 
magnetic moment of the muon, {g — 2)^, from various energy regimes are determined and, 
with the addition of experimental data, final values for a(M|) and = {g — 2)^/2 are 
determined. 



2 Running of the QED Fine Structure Constant 

The running of the electromagnetic fine structure constant a{s) is governed by the renor- 
malized vacuum polarization function, 11^ (s). For the spin 1 photon, n^(s) is given 
by the Fourier transform of the time-ordered product of the electromagnetic currents 
i^^{s) in the vacuum (gV - ^V'') n^(9^) = i S d'^x e''^'' {Q\T{i^^{x)il^{<d))\Q) . With 
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Aa{s) — 47raRe [n.y(s) — 11^(0)], one has 

a{s) = 



a(0) 



Aais] 



(2) 



where 47ra(0) is the square of the electron charge in the long- wavelength Thomson limit. 
The contribution Aa{s) can naturally be subdivided in a leptonic and a hadronic part. 



The leading order leptonic contribution is given by 



on 



VTLc 



Ml, 



314.2 X 10" 



(3) 



Using analyticity and unitarity, the dispersion integral for the contribution from the light 
quark hadronic vacuum polarization Aahad(Afz) reads 



Aahad(M|) 



Ml 



47?^ a 



Re / ds 



0-had(S 



Ml 



(4) 



where a"had(s) = IGvr a ■ Imn^(s) is given by the optical theorem, and Imll^ stands 

for the absorptive part of the hadronic vacuum polarization correlator. Through Eq. (|l|), 
the above dispersion relation can be expressed as a function of -R(s): 



A«had(M|) 



Re / ds 



R{s) 



Ami 



s{s - Ml) - ie 



(5) 



Employing the identity l/{x' — x — ie)^^o = P{l/{x' — x)} + m5{x' — x), the integrals (^) 
and (13) are evaluated using the principle value integration technique. 



3 Muon Magnetic Anomaly 

It is convenient to separate the prediction from the Standard Model into its different 
contributions 

= af ° + a^;' + a^;^^ , (6) 

where a|^™ = (11 658 470.6 ± 0.2) x 10^^° is the pure electromagnetic contribution (see [O 
and references therein), a^"^ is the contribution from hadronic vacuum polarization, and 
aj^^^"^ = (15.1 ± 0.4) X 10~^° |Tl], |1^, |T^ accounts for corrections due to the exchange of 



the weak interacting bosons up to two loops. 

Equivalently to Aahad(M|), by virtue of the analyticity of the vacuum polarization 
correlator, the contribution of the hadronic vacuum polarization to can be calculated 



via the dispersion integral |14 



Ami 
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Here K{s) denotes the QED kernel 



K{s) = (^1 - y j + (1 + (^1 + ^) (^ln(l + x) - X + y j + Inx (8) 

with a; = (1 — + /S^^) and /5 = (1 — 4m^/s)^/^. The function -ft'(s) decreases 

monotonically with increasing s. It gives a strong weight to the low energy part of the 
integral @. About 91% of the total contribution to a^^'^ is accumulated at cm. energies 
^/s below 2.1 GeV while 72% of a^^'^ is covered by the two-pion final state which is 
dominated by the p(770) resonance. Data from vector hadronic r decays published by the 
ALEPH Collaboration provide a very precise spectrum of the two-pion final state as well 
as new input for the lesser known four-pion final states. This new information improves 
significantly the precision of the a^^*^ determination [||]. 



Theoretical Prediction of R(s 



The optical theorem relates the total hadronic width at a given energy-squared Sq to the 
absorptive part of the photon vacuum polarization correlator 

R{sq) = 127rImn(so + «e) . (9) 

Perturbative QCD predictions up to next-to-next-to leading order are available for the 
Adler D-function [|16| which is the logarithmic derivative of the correlator 11, carrying all 
physical information: 

D{s) = -Utt's^^ . (10) 
as 



This yields the relation 



\s\=so 



where the contour integral runs counter-clockwise around the circle from s = Sq — ie to 
s = Sq + ie. Choosing the renormalization scale to be the physical scale s, additional 
logarithms in the perturbative expansion of D are absorbed into the running coupling 
constant ^^(s). The (massless) NNLO perturbative prediction of D reads then |jl^ 
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Dp{-s) = NcY.Q^ 
f 




+ O {at{s)) 



(12) 



where Nc = 3 for SU (3)c and Qf is the charge of the quark /. The coefficients are do = 1, 
di = 1.9857-0.1153 ri/, rfs = rf2+/5ovrV48, /3o = ll-2n//3 and = -6.6368-1.2001 ri/- 
0.0052 - 1.2395 {T.fQff/NcY.fQ) with nf the number of involved quark flavours^. 
The running of the strong coupling constant as{s) is governed by the renormalization 

^ The negative energy-squared in _Dp(— s) of Eq. ( p^ is introduced when continuing the Adler function 
from the spacehke EucUdean space, where it was originally defined, to the timelike Minkowski space by 
virtue of its analyticity property. 
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group equation (RGE), known precisely to four-loop level ||TB[. 



Using the above formalism, R{so) is easily obtained by evaluating numerically the 
contour-integral (0). The solution is called contour- improved fixed-order perturbation 
theory (FOPTci) in the following. Another approach, usually chosen, is to expand as{s) 
in Eq. (0) in powers of as{so) with coefficients that are polynomials in ln(s/so): 



a.As 



TT 



-/3oln— 

71 4: So 




8 So 
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T^/52ln---^/3oAln^- + -^/3o^ln3- 
128 So 64 So 64 sq 




(13) 



Inserting the above series with the Dp-function in Eq ([TT|) and keeping terms up to oder 
leads to the expression 



/2(so) 



f 



1+do^+d, 



TT 



n 



+ d2 



7C 




O (atiso)) 



(14) 



where the difference to Dp is of order only. The solution (p^ ) will be referred to as 
fixed-order perturbation theory (FOPT). There is an intrinsic ambiguity between FOPT 



and FOPTci- The numerical solution of the contour-integral (11) involves the complete 
(known) RGE and provides thus a resummation of all known higher order logarithmic 
terms of the expansion (|T3|) (see Ref. for comparison). Unfortunately, it is unclear if 
the resummation does not give rise to a bias of the final result. 



Quark Mass Corrections 

Quark mass corrections in leading order are suppressed as ~ m^(s)/s, i.e., they are small 
sufficiently far away from threshold. Complete formulae for the perturbative prediction 




containing quark masses are provided in Refs. ||2y, ^ for the correlator n(s) and R{s) 
up to order exactly and numerically using Pade approximants to order a^. We will 
learn from the numerical analysis that it suffices for the required level of accuracy to use 
the following expansion as an additive correction to the Adler D-function [^] 

D^U-s) = -NcY.Q''^(^Q + '^8^ + (294.8 -12.3nf)[^] | . (15) 

The running of the quark mass mf{s) is obtained from the renormalization group and is 
known to four-loop level P^ . 

When using perturbative QCD at low energy scales one has to worry whether contri- 
butions from nonperturbative QCD could give rise to large corrections. The break down 
of asymptotic freedom is signalled by the emergence of power corrections due to nonper- 
turbative effects in the QCD vacuum. These are introduced via non-vanishing vacuum 
expectation values originating from quark and gluon condensation. It is convenient to use 
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the Operator Product Expansion (OPE) |2^, in low energy regions (or near quark 
thresholds), where nonperturbative effects come into play. One thus defines 



D{-s) 



E 

D=0,2,4,, 



(16) 



dimO=D 



where the Wilson coefficients C{s,^) include short- distance effects and the operator 
{Ot>{^)) collect the long-distance, nonperturbative dynamical information at the arbi- 
trary separation scale fi. The operator of Dimension D = {D = 2) is the perturbative 
prediction (with mass correction). The D = A operator are linked to the gluon and quark 
condensates. Effective approaches together with the vacuum saturation assumption are 
used to compact the large number of dynamical D = 6 {D = 8) operator into one phe- 
nomenological operator {Oq) {{Os}). The nonperturbative addition to the Adler function 
reads then 



■5-^ of 27r^ / 11 a,(s 





^GG 



18 TT / 



27 TT 



E 



{rrikqkqk) 



(17) 



with the gluon condensate, {{as/7!')GG), and the quark condensates, {nifqfqf). The latter 
obey approximately the PCAC relations 



(m„ + nid) {uu + dd) 
ms{ss) 



p2 2 



where = (92.4 ± 0.26) MeV |2^ is the pion decay constant. The complete dimension 
D = 6 and D = 8 operator are parametrized phenomenologically using the vacuum ex- 
pectation values (Oe) and (Og), respectively. Note that in zeroth order a^, i.e., neglecting 
running quark masses, nonperturbative dimensions do not contribute to the integral in 
Eq. ([TI|) . Thus in the formula presented in Eq. ([T8|) only the gluon and quark condensates 
contribute to R via the logarithmic s-dependence of the terms in first order a^. 

The total Adler D-function then reads as the sum of perturbative, mass and nonper- 
turbative contributions: 



D(s) = Dp(s) + 



^NP(S) 



(19) 



Uncertainties of the Theoretical Prediction 

Looking at Eq. (|T^) it is instructive to subdivide the discussion of theoretical uncertainties 
into three classes: 
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[i) The perturbative prediction. The estimation of theoretical errors of the perturbative 
series is strongly linked to its truncation at finite order in as- Due to the incomplete 
resummation of higher order terms, a non- vanishing dependence on the choice of the 
renormalization scheme (RS) and the renormalization scale is left. Furthermore, one 
has to worry whether the missing four-loop order contribution ds{as/ vr)^ gives rise to 
large corrections to the series (0). On the other hand, these are problems to which 
any measurement of the strong coupling constant is confronted with, while their 
impact decreases with increasing energy scale. The error on the input parameter 
itself therefore reflects to some extent the theoretical uncertainty of the perturbative 
expansion in powers of as- 

Let us use the following, intrinsically different as determinations to benchmark 
our choice of its value and uncertainty. A very robust as measurement is obtained 
from the global electroweak fit performed at the Z-boson mass where uncertainties 
from perturbative QCD are rather small. The value found is q;s(M|) = 0.120 ± 
0.003 0. A second precise as measurement is obtained from the fit of the OPE 
to the hadronic width of the r and to spectral moments 0. The measurement is 
dominated by theoretical uncertainties, from perturbative origin. In order to ensure 
the reliability of the result, i. e., the applicability of QCD at the r mass scale, spectral 
moments were fitted analogously to the analysis presented in the following section. 
The nonperturbative contribution was found to be lower than 1%. Additional tests 
in which the mass scale was reduced down to 1 GeV proved the excellent stability 
of the as determination. The value reported by the ALEPH Collaboration [^j is 
as(M|) = 0.1202 ± 0.0026. A third, again different approach is employed when 



using lattice calculations fixed at bb states to adjust as- The value given in Ref. 
is a,(M|) = 0.117 ±0.003. 

The consistency of the above values using quite different approaches at various 
mass scales is remarkable and supports QCD as the theory of strong interactions. 
To be conservative, we choose q;s(M|) = 0.1200 ± 0.0045 as central value for the 
evaluation of the perturbative contribution to the Adler D-function. 

Even if it is in principle contained in the uncertainty of Aa<j(M|) = 0.0045, we 
furthermore add the total difference between the results obtained using FOPTci 
and those from FOPT as systematic error. 

(a) The quark mass correction. Since a theoretical evaluation of the integral (^) is only 
applied far from quark thresholds, quark mass corrections -Dmass are small. Without 
loss of precision, we take half of the total correction as systematic uncertainty, i.e., 
we add Dmass ± -Dmass/2 iu Eq. ([T^). 

(Hi) The nonperturbative contribution. In order to detach the measurement from theoret- 
ical constraints on the nonperturbative parameters of the OPE, we fit the dominant 
dimension D = 4,6,8 terms by means of weighted integrals over the total e~^e~ low 
energy cross section. Again, without loss of precision, we take the whole nonpertur- 
bative correction as systematic uncertainty, i.e., we add D^p ± D^p in Eq. ([T9|). 

Another sources of tiny uncertainties included are the errors on the Z-boson and the top 
quark masses. 
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5 Spectral Moments 



Constraints on the nonperturbative contributions to R{s) from theory alone are scarce. 
It is therefore advisable to benefit from the information provided by the explicit shape 
of the hadronic width as a function of s in order to determine the magnitude of the OPE 
power terms at low energy. We consequently define the following spectral moments 



^ y ds f s 



R'\so)= / - 1-- (- Ris) , (20) 

J So V Sq/ \Sq/ 

where the factor (I — s/sq)^ squeezes the integrand at the crossing of the positive real axis 
where the validity of the OPE is questioned. Its counterpart (s/so)' projects on higher 
energies. The new spectral information is used to fit simultaneously the phenomenological 
operators {{as/TT)GG), {Oq) and (Cg), a procedure which requires at least 4 — better 5 — 
input variables considering the intrinsic strong correlations between the moments which 
are reinforced by the experimental correlations and by the correlations from theoretical 
uncertainties. 

To predict theoretically the moments, one uses the virtue of Cauchy's theorem and 
the analyticity of the correlator n(s), since a direct evaluation of the integral (^) in the 
framework of perturbation theory (and even OPE) is not possible. With the relation (^, 
Eq. ( pO|) becomes 

R'\so) =6m i -{l- -t (- j n(.) , (21) 
J So V So/ Vso/ 

\s\=so 

and with the definition ( p!0D of the Adler D-function one further obtains after integration 

by parts 



1 r ds 
\ f 

s|=so 

where -D(s) is obtained from Eq. (|19|). 



sl=sn L 



D{s) , (22) 



6 Data analysis and Determination of the Moments 

Due to the suppression of nonperturbative contributions in powers of the energy scale s, 
the critical domain where nonperturbative effects may give residual contributions to R{s) 
is the low-energy regime with three active flavours. We thus choose the energy scale So 
of the fit equal to the energy scale, where the theoretical evaluation of R{s) shall start. 
As demonstrated in isovector vector r decays 0, the scale ^s = Mr is an appropriate 
scale where nonperturbative effects are present, but essentially controlled by the OPE. 
In our case we manipulate isovector and isoscalar vector hadronic final states, i.e., more 
inclusive data, and might expect smaller nonperturbative contributions. We therefore set 
the energy cut to ^sq = 1.8 GeV. Up to this energy, R[s) is obtained from the sum of 
the hadronic cross sections exclusively measured in the occurring final states. 



7 



Moments {k,l) 


Data 


A(stat. + sys.) 


A(theo.) 


(2,0) 


0.770 


0.013 


0.024 


(2,1) 


0.204 


0.005 


0.012 


(2,2) 


0.073 


0.003 


0.001 


(2,3) 


0.035 


0.002 


0.001 


(2,4) 


0.020 


0.001 


0.001 




Moments [k, I) 


(1,0) 


(1,1) (1,2) 


(1,3) (1,4) 


(1,0) 


1 


0.98 0.75 


0.60 0.49 


(1,1) 




1 0.79 


0.67 0.57 


(1,2) 




1 


0.97 0.91 


(1,3) 






1 0.99 


(1,4) 






1 



Table 1: Spectral moments measured at ^,/sQ = 1.8 GeV with experimental and theoreti- 
cal errors. Below, the corresponding correlation matrix containing the quadratic sum of 
experimental and theoretical covariances. 

The data analysis follows exactly the line of Ref. In addition to the e~^e~ annihi- 
lation data we use spectral functions from r decays into two- and four final state pions 
measured by the ALEPH Collaboration p9|. Extensive studies have been performed in 
Ref. in order to bound unmeasured modes, such as some KK or the 7r+7r~47r° final 
states, via isospin constraints. We bring attention to the straightforward and statisti- 
cally well-defined averaging procedure and error propagation used in this paper as in the 
preceeding one, which takes into account full systematic correlations between the cross 
section measurements. All technical details concerning the data analysis and the integra- 
tion method used are found in Ref. 

The experimental determination of the spectral moments ( pOD is performed as the sum 
over the respective moments of all exclusively measured e~^e~ final states (completed by r 
data). We chose the moments = 2, Z = 0, ■ ■ ■ , 4, in order to collect sufficient information 
to fit the three nonperturbative degrees of freedom. Neglecting the s-dependence of the 
Wilson coefficients in Eq. (|16D, the respective nonperturbative power terms contribute to 
the following moments: the dimension D = 4 term contributes to / = 0, 1, the D = 6 term 
to / = 0, 1, 2 and the D = 8 term contributes to the / = 0, 1, 2, 3 moments. The / = 4 
moment receives no direct contribution from any of the considered power terms. However 
its use is not obsolete, since it constrains the power terms through its correlations to the 
other moments. Table |I| shows the measured moments together with their (statistical 
and systematic) experimental and theoretical errors. Additionally given is the sum of the 
experimental and theoretical correlation matrix as it is used in the fit. 

Using as input parameters as(M|) = 0.1200 ± 0.0045, yielding for three flavours 
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Ag^ = (372 ±76) MeV), and for the mass of the strange quarkgm,(l GeV) = 0.230 GeV, 
while setting = = 0, the adjusted nonperturbative parameters are 

(0.037 ±0.019) GeV^ , 
-(0.002 ±0.003) GeV^ , 

(0.002 ± 0.003) GeV^ , (23) 

with x^/d.o.f = 1.5/2. The correlation coefficients between the fitted parameters are 
p{{{ajTi)GG),{0,)) = -0.41, p(((«,/7r)GG),(08)) = 0.55 and p{{0,),{0^)) = -0.98. 
As a test of stability we have additionally fitted the A; = 2, / = 0, . . . , 4 spectral moments 
at a/so = 2.1. The difference between the results of this fit and Eq. (pSf ) is included 
into the parameter errors given in Eq. (|23|). With these values, the corrections to R{sq) 
at a/sq = 1.8 GeV amount to 0.16% from the strange quark mass and < 0.1% from the 
nonperturbative power terms. The gluon condensate can be compared to the standard 
value obtained from charmonium sum rules, {{as/n)GG) = (0.017 ± 0.004) GeV^ P3| , 
which lies below our value. However, another estimation using finite energy sum 
rule techniques on e~^e~ data gives the value of {{as/i^)GG) = (0.044 1o;q2i) GeV^ in 
agreement with the result (^). Fitting the moments when fixing the dimension D = 
6, 8 contributions reduces the gluon condensate to 0.010 ± 0.002. One may additionally 
compare the fitted dimension D = 6,8 operator to the results obtained from the r vector 
spectral functions, keeping in mind that only the isovector amplitude contributes in this 
case and thus the more inclusive isoscalar plus isovector moments from e~^e~ annihilation 
are expected to receive smaller nonperturbative contributions^. Reexpressing the results 
of Ref. P in terms of the definition adopted in Eq. ( [T7| ) we obtain {Oq)i=i = — (0.0042 ± 
0.0006) GeV^ and (08)/=i = (0.0062 ± 0.0007) GeVl 

As a cross check of the spectral moment analysis, we fit the three nonperturbative 
parameters and a^. The theoretical error applied reduces essentially to the theoretical 
uncertainties of the QCD perturbative series estimated in Ref. to be Aas(M|) = 0.0023 
at the T mass scale which is of the same magnitude as the scale ^/sq = 1.8 GeV used here. 
The fit of the A; = 2, / = 0, . . . , 4 moments yields a,(M|) = 0.1205 ± 0.0053 (x^/d.ol = 
1.4/1), in agreement with the values from other analyses cited above. The values for 
the gluon condensate and the higher dimension contributions are consistent with those 
of Eq (^). Fitting ^^(so) and the gluon condensate at ^/so = 2.1 GeV when fixing the 
dimension D = 6,8 contributions at the values of Eq. (p3|) yields as(M|) = 0.121 ± 0.011 

^ The mass of the strange quark used in this analysis takes one's bearings from the recent experimental 
determination using the hadronic width of r decays into strange final states, Rt.s, reported by the 
ALEPH collaboration [|§ to be ms(l GeV) = 235t^^ MeV. The error on this mass has no influence on 
the present analysis since, conservatively, 50% of the total mass contribution given in Eq. ( p^ is taken 
as corresponding systematic uncertainty. 

^ The vacuum saturation hypothesis offers a relation between the dimension D — 6 contribution and 
the light quark condensates Using the formulae ( p7| ) and (|l^) one has 

004 
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GG 

(Os) 
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and {{as/ 11)00) = (0.042 ± 0.003) GeV. The consistency of these results with those 
obtained at ^/so = 1.8 GeV supports the stabihty of the OPE approach within the energy 
regime where it is apphed. Varying the cm. energy sq and the weights k, I of the moments 
used to fit the nonperturbative contributions gives a measure of the compatibihty between 
the data and the OPE approach. Differences found between the fitted parameters are 
included as systematic uncertainties in the errors of the values (PBI). 



7 Evaluation of Aahad(M|) and aj;^^ 

In order to get the most reliable central value of the perturbative R{s) prediction which 
enters the integrals (H) and (|^), we use the whole set of formulae given in Ref. includ- 



ing mass corrections up to order a^. Despite this theoretical precision, the uncertainties 
keep conservatively estimated as described in Section |^. 

Since deeply nonperturbative phenomena are not predictable within the OPE ap- 
proach, we use experimental data to cover energy regions near quark thresholds. The 
low energy results for Aahad(Af|) and a^^'^ from Ref. [0] including r data are taken for 
y/s < 1.8 GeV. The narrow uj, 0, J/ip and T resonances are parametrized using relativis- 
tic Breit-Wigner formulae as described in Refs. p2| , 0]. In addition, R measurements of 



the continuum contributions in the environment of the cc threshold are taken from the 
experiments specified in Ref. [0] . The technical aspects of the integration over data points 
are also discussed in Ref. ||l|. No continuum data are available at bb threshold energies in 
the range of 10.6 GeV < E < 12 GeV. A recent analysis of the T system |^ however 



showed that essentially within the perturbative approach of Ref. ||2T| it is possible to pre- 
dict weighted integrals over the bb states. The uncertainty of this approach corresponds 
to an estimated error of Aas(M|) = 0.008. We therefore use this uncertainty of for 
the QCD prediction at bb threshold energies. Included is a small uncertainty originating 
from scale ambiguities when matching the effective theories of four and five flavours. 

For the theoretical evaluation of the integrals (^), (0) (via Eq. ([TT|)), we use the fol- 
lowing variable settings: 

Mz = (91.1867 ±0.0020) GeV 0, 
a(0) = 1/137.036 , 
as{Mi) = 0.1200 ± 0.0045 , 

mu = rrid = , 

m,(l GeV) = 0.230 GeV , 

Tnc{mc) = 1.3 GeV , 

fnb{^b) = 4.1 GeV , 

mt{mt) = (175.6 ± 5.5) GeV , 

and the values (^) for the nonperturbative contributions. There are no errors assigned 
to the light quark masses, since the half of the total quark mass correction is taken as 
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-E'min - -E'max (GeV) 


Aahad( 
Data 


M|) X 10^ 

Theory 


cr 


^had X ^QlO 


Ami - 1-8 


56.9 ± l.l(*) 






636.49 ± 7.41 (D) 


1.8 - 3.700 


32.4 ±3.1 


24.50 ± 0.33(*) 


2.5 


33.84 ±0.53 (T) 


^^(3770) 
3.700 - 5.000 


0.29 ±0.08^ 
15.8 ± I.7W 


17.06 ±0.58 


0.5 


0.17 ±0.02 (D) 
6.93 ±0.62 (D) 


5.000 - 10.500 


39.9 ± 1.4 


41.55 ±0.40(*) 


1.0 


7.43 ± 0.09 (T) 


T(4S, 10860, 11020) 
10.500 - 12.000 


0.38 ± 0.10 
7.6 ±0.5 


8.19±0.32W 


0.4 


0.55 ±0.03 (T) 


12.000 - 40.000 


75.2 ±2.7 


77.96 ± 0.30(*) 


1.0 


1.64 ±0.02 (T) 


40.000 - oo 




41.98 ±0.22W 




0.16 ±0.00 (T) 


J/^(1S,2S) 


9.68±0.68W 






7.80 ± 0.46 (D) 


T(1S,2S,3S) 


0.98 ±0.15^ 






0.09 ±0.01 (D) 


4m^ - oo 


277.8 ± 2.2exp ± 1.4theo 


695.1 ± 7.5exp ± 0.7theo 



Value used for final result of AahadiM^) (last line). 



Table 2: Contributions to Aa;had(^z) ^'^'^ ^'^ ^^^'^ from the different energy regions. 
The "a" column gives the standard deviation between the experimental and theoretical 
evaluations of Aahad(^z)- in the last column stands for evaluation using 

Data/Theory. 

systematic uncertainty. The error on nit is needed in order to estimate the systematic un- 
certainty of the matching scale when turning from five to six flavours. Table |^ shows the 
experimental and theoretical evaluations of Aahad(^z) ^^^'^ respective energy 

regimes. The upper star denotes the values used for the final summation of Aahad(^z) 
given in the last line. At ~ 3.8 GeV a non-resonant DD production might contribute to 
the continuum^. We therefore use experimental data to cover energies from 3.7-5.0 GeV. 

A 20% correlation is assumed between the analytic evaluations of the narrow reso- 
nances where in each case a Breit-Wigner formula is applied. The theoretical errors are by 
far dominated by uncertainties from as and the difference FOPT/FOPTci. For instance, 
the first energy interval where theory is applied {E G {1.8 — 3.7 GeV}) receives the error 



^ Such a contribution must be tiny since it is suppressed by its form factor and the (1 — 4A/^/s)^/^ 
threshold behaviour of a pair of spin particles. 
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contributions 



Aahad(M|) 



X 10^ = 24.502 ± 0.332 



0.243 - Aos 

0.223 — perturbative prediction 

0.036 — quark mass correction 

0.003 — nonperturbative parts 

< 0.001 - AMz 



where the small contributions from quark masses and nonperturbative dimensions show 
that the perturbative QCD calculation is very solid here. Remember that only Ins- 
dependent nonperturbative terms contribute to R{s). Theoretical errors of different en- 
ergy regions are added linearily except the uncertainties at bb threshold that are (partly) 
from individual origin so that a 50% correlation to other energy regimes is estimated here. 
Looking at Table |^ one notices the remarkable agreement between experimental data and 
theoretical predictions of Aahad(^z) even in the cc quark threshold regions where strong 
oscillations occur. The experimental results of R{s) and the theoretical prediction are 
shown in Fig. |l|. The shaded bands depict the regions where data are used instead of 
theory to evaluate the respective integrals. Good agreement between data and QCD is 
found above 8 GeV, while at lower energies systematic deviations are observed. The R 



measurements in this region are essentially provided by the 772 |3^ and MARK I 
collaborations. MARK I data above 5 GeV lie systematically above the measurements of 
the Crystal Ball 



and MDl [pi] Collaborations as well as the QCD prediction. 



8 Results 



According to Table the combination of the theoretical and experimental evaluations of 
the integrals (H) and (0) yield the final results 



A«had (M|) 


= (277.8 ± 2.2exp ± 1.4theo) X 10-4 


a-\Mi) 


= 128.923 ± 0.030exp ± 0.019theo 


had 


= (695.1 ± 7.5exp ± 0.7thco) X 10-^° 




= (11 659 164.6 ± 7.5exp ± 4.1theo) x lO^i" 



(24) 



The total value includes an additional contribution from non-leading order hadronic 
vacuum polarization summarized in Refs. |3^, |I| to be a^^'^^a / n)^] = (—16.2 ± 4.0) x 
10"^*^, where the error originates essentially from the uncertainty on the theoretical evalu- 
ation of the light-by-light scattering type of diagrams ^ . Fig. ^ shows a compilation 
of published results for the hadronic contribution Aahad(^z)- Some authors give the 
hadronic contribution for the five light quarks only and add the top quark contribution 
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0^ 




e e'^hadrons 
— QCD 

J I \ I I \ I I I L 



_l I I \ I I I \ I I I \ I I I I I L 



* MARK I * MDl 

A PLUTO * JADE 

A LENA o MARK J 

■ Crystal B. □ CESR 

I I \ I I I \ I I I \ L_ 



7 



10 11 12 13 14 

Vs (GeV) 



Figure 1: Inclusive hadronic cross section ratio in e'^e' annihilation versus the cm. 
energy ^/s. Additionally shown is the QCD prediction of the continuum contribution as 
explained in the text. The shaded areas depict regions were experimental data are used 
for the evaluation o/ Acthadl^l) ^'^^ '^^'^ addition to the measured narrow resonance 
parameters. The exclusive e~^e~ cross section measurements at low cm. energies are 
taken from DM1,DM2,M2N,M3N,0LYA,CMD,ND and t data from ALEPH (see Ref. ^ 
for detailed information) . 
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Figure 2: Comparison of estimates of and Aa^adiM^] 

Refs. m ii,ii,ii,ii,n. 



The values are taken from 



of AQ;top(M|) = —0.6 X 10~^ separately. This has been corrected for in the figure. The 
present inclusion of theoretical evaluations yields an improvement in the Aahad(Af|) pre- 
cision of a factor of two. The improvement of the precision on a^^'^ with respective to the 
analysis presented in Ref. amounts to 20%. 

The analysis using the spectral moments and also the measurement of at the r mass 
scale showed concordantly that the OPE can be safely applied in order to predict inte- 
grals over inclusive spectra at relatively low energy scales. Nonperturbative contributions 
are indeed tiny and well below the envisaged precision. This approach leads to a large 
improvement in precision, as compared to the method used in Ref. |^3| where perturbative 
QCD was assumed to be valid only above 3 GeV. Also the value used for as(M|) in the 
latter analysis was less precise than the current value. Finally, better and more complete 
experimental information, in particular at low energy, is used in the present work. Exper- 
imental data on R are still necessary at very low energies and at the cc quark production 
threshold. Uncontrolled nonperturbative effects spoil the hadronic spectra and make the 
OPE approach unreliable. Thus, the final results on both, Aahad(^z) ^^"^ (^^^'^^ ^ire still 
dominated by experimental uncertainties, where in particular the unmeasured tt^ti^Ati^ 
and KKtttt low energy final states, which must be conservatively bound using isospin 
symmetry as well as the two-pion final state in the latter case, contribute with large 
errors. Also more precise data of the cc continuum at threshold energies are needed. 

We use the improved precision on a(M|) to repeat the global electroweak fit in or- 
der to adjust the mass of the Standard Model Higgs boson, Mniggs- As electroweak and 
heavy flavour input parameters we use values measured by the LEP, SLD, CDF, DO, 
CDHS, CHARM and CFFR collaborations, which have been collected and averaged in 
Ref. . The prediction of the Standard Model is obtained from the ZFITTER electroweak 
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library The fit adjusts Mz, Mtop and a{Mz) which are allowed to vary within their er- 
rors. Freely varying parameters are and Mniggs- We obtain as{M^) = 0.1198 ± 0.0031 
in agreement with the experimental value of 0.122 ± 0.006 from the analyses of QCD 



observables in hadronic Z decays at LEP. The fitted Higgs boson mass is 129lg2^ GeV 
with = 16.4/15, compared to 105lg2^ GeV when using the previous value of a(M|) 
from Ref. An additional error of 50 GeV should be added to account for theoretical 



uncertainties Doing so, we obtain an upper limit for Mniggs of 398 GeV at 95% CL. 

Fig. 1^ depicts the variation of Ax^ as a function of the Higgs boson mass for the new 
and previously used values of a(M|) [|l|. 

Table 1^ shows the dominant uncertainties of the input values of the Standard Model fit 



expressed in terms of Asin^^J^ . This work reduces the uncertainty of a(M|) on Mnig, 
below the uncertainties from the experimental value of 8110^9]^^^ or from the Standard 
Model, i.e., theoretical origin. 



9 Conclusions 

We have reevaluated the hadronic vacuum polarization contribution to the running of the 
QED fine structure constant, a{s), at s = M| and to the anomalous magnetic moment of 
the muon, a^. We employed perturbative and nonperturbative QCD in the framework of 
the Operator Product Expansion in order to extend the energy regime where theoretical 
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Experiment 
a(M|) 

rrit 
Theory 

^Higgs 



0.00023 
0.00023 ^^'^'^ 0.00010 
0.00018 
0.00014 

Fit 



0.00160 [Mmggs = 60-1000 GeV] ^ Mniggs 



I29tlf GeV 



Table 3: Dominant uncertainties of input values of the Standard Model electroweak fit 
expressed in terms o/ A sin^^^'^'* 0, |^. 

predictions are reliable. In addition to the theory, we used data from e^e~ annihilation and 
r decays to cover low energies and quark thresholds. The extended theoretical approach 
reduces the uncertainty on Aahad(^z) ^y more than a factor of two. Our results are 
Aahad(M|) = (277.8 ± 2.6) x 10~\ propagating a'\0) to a-\M^) = (128.923 ± 0.036), 
and a^'^'^ = (695.1 ± 7.5) x 10^° which yields the Standard Model prediction = 
(11659164.6 ± 8.5) x 10~^°. The new value for q;(M|) improves the constraint on the 
mass of the Standard Model Higgs boson to Mniggs = 129lg2^ GeV. 
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